Effective super Tonks-Girardeau gases as ground states of strongly attractive 
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In the strong interaction limit, attractive fermions with iV-component hyperfine states in a one- 
dimensional waveguide form unbreakable bound cluster states. We demonstrate that the ground 
state of strongly attractive SU(iV) Fermi gases can be effectively described by a super Tonks- 
Girardeau gas-like state composed of bosonic cluster states with strongly attractive cluster-cluster 
interaction for even N, and a Fermi duality of a super Tonks-Girardeau gas-like state composed of 
fermionic cluster states with weakly interacting cluster-cluster p-wave interaction for odd N. 

PACS numbers: 03.75.Ss, 05.30.Fk 
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I. INTRODUCTION 



The experimental progress on trapping ultracold atoms 
in tightly confined waveguides in a well-controlled way 
[l]-[4| has stimulated intensive study of the physical prop- 
erties of one-dimensional (ID) quantum gases. The 
effective interaction strength between atoms in a ID 
waveguide can be tuned via Feshbach resonance or 
confinement- induced resonance pj, leading to the exper- 
imental realization of Tonks-Girardeau (TG) gases [3, La]. 
The TG gas describes the strongly repulsive Bose gas 
0, 0|- Starting from the TG gas and then switching 
the interaction between atoms from strongly repulsive 
to strongly attractive, the experimental realization of a 
ID super Tonks-Girardeau (STG) gas of bosonic Cesium 
atoms was reported very recently |8(. In contrast to the 
TG gas, the STG gas describes a gas-like phase of the 
attractive Bose gas f9l-fl"lj|. which is metastable against 
falling into its cluster-type ground state [TH , despite the 
fact that the interaction between atoms is strongly at- 
tractive (HQ. 

Although the STG gas realized in Q is a metastable 
highly excited state of the attractive Bose gas, in recent 
theoretical work it was found that the ground state of 
a strongly attractive spin- 1/2 Fermi gas can be effectively 
described by the STG gas. Intuitively, two fermions with 
opposite spins form a tightly bound state and the bound 
pairs of fermions can be viewed as composite bosons with 
a mass of 2m f ■ The effective interaction between the 
composite bosons is also attractive with the interaction 
strength given by cb — 2cp [HJ EH- Conversely, the 
ground state of the bound Fermi pairs is described by 
the STG phase of attractive bosons [l5|. Very recently, 
multi-component Fermi gases have attracted considerable 
interest |17H19| due to the novel existence of different 
sizes of molecules. Here we consider the interesting ques- 
tion of whether the ground state of a strongly attractive 
multi-component fermionic system can be also effectively 
described by a STG gas of multi-particle bound states. 



A positive answer is confirmed by explicit identification 
of the general mapping relation between the attractive 
SU(iV) Fermi gas and the STG gas. 



II. ATTRACTIVE iV-BOUND FERMIONS 

We consider a delta-function interacting system com- 
posed of Np atomic attractive fermions with equal mass 
mp which occupy N hyperfine levels with identical par- 
ticle numbers N l = N N = Np /N (i = l,...,N) and 
constrained by periodic boundary conditions to a line of 
length L. If the interactions are spin independent, the 
Hamiltonian reads 



d 2 



Imp dx 2 



9F 



^5{xi-x 3 ), (1) 



where gp = — 2h 2 /(mp af D ) is the interaction strength 
with af D the ID effective s-wave scattering length. Al- 
though the interactions in Hamiltonian (|TJ) are repre- 
sented in terms of delta interaction, the exchange anti- 
symmetric wavefunction for fermions gives the restriction 
that interactions among the same fermion level are pro- 
hibited. Different symmetries of the wavefunction pro- 
duce different Bethe ansatz equations even if the Hamil- 
tonian has the same form as Eq. ([TJ [2(| l2lj |- For 
simplicity, we use the dimensionless coupling constant 
If = cp /np with density np — Np/L and cp — —2/af D . 
In the strongly attractive limit (gp — > — oo), atoms form 
tightly bound states with each bound state composed of 
N fermions in different hyperfine states [22| , No tightly 
bound state with more than N fermions can be formed 
in an V-component Fermi gas due to the Pauli exclusion 
principle. 

If N is even, the tightly bound state can be viewed 
as a composite boson with a mass mg = Nmp. In this 
work, we find that the ground state of the strongly at- 
tractive SU(Af) Fermi gas with TV even can be effectively 
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described by a STG gas of attractive composite bosons. 
This can be viewed as a direct generalization of the SU(2) 
Fermi gas result On the other hand, if N is odd, the 
tightly bound state can only be viewed as a composite 
fermion. Obviously, strongly attractive SU(AT) fermions 
with N odd cannot be mapped to an effective bosonic 
gas as for the two-component case [Hj]. Nevertheless, we 
shall show that it can be mapped to a spinless Fermi 
gas with weakly interacting p-wave interaction. Due to 
a general Fermi-Bose mapping [23| , eigenstates of a spin- 
less Fermi gas with p-wave interaction of any strength 
can be mapped to those of a ID Bose gas with delta- 
function interactions. Therefore the gas-like state of a 
weakly repulsive p-wave Fermi gas can be viewed as the 
Fermi correspondence of the STG phase of a strongly 
attractive Bose gas. 

Before construction of the mapping between the 
ground state of the strongly attractive SU(A r ) Fermi gas 
and the STG phase of composite bosons or fermions, we 
first discuss the solution of the SU(AT) Fermi gas which 
is exactly solvable by the Bethe ansatz method. The 
eigenvalues of Hamiltonian (jTJ) are given by 



E = 



.2 

2m F ^ '" j 



(2) 



with quasi-momentum kj determined by the Bethe 

ansatz equations (BAEs) [20W22Tj 



exp(ifcjL) = Yl 



r kj A a ^ + ICp 
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for j = 1,...,N F , a = 1, 
1. We have denoted M = Np and A 



(0) 



(4) 

,N- 
The 



parameters are the spin rapidities. The quantum 

numbers are given by Mi = (N — I) Nn and c' F = cf/2. 

For strongly attractive attraction, i.e., L\cp\ 3> 1, the 
BAEs permit different sizes of charge bound state. As 
a consequence of the Pauli exclusion principle and the 
SU(N) symmetry, there is no tightly bound state with 
more than N fermions for the SU (AT) Fermi gas [lij . For 
the ground state, there are equal numbers of particles in 
each hyperfme spin state. In this state, the charge bound 
state in fc-space is of the form 

k qd = Af-D + (N + 1 - 2j) c' F + (i6\c F \) , (5) 

for j = 1,2,..., N and q = 1,2, N F /N. The spin 
rapidities form a certain pattern of spin string solutions. 
For the ground state, each charge bound state kgj with 

different q is accompanied by a spin string |a^1| of 



length N — lwith a = 1,2, ... , N—l, a spin string jA^a j 
of length N — 2 with a = 1,2, . . . , N — 2, and so on, until 
the real root A[ N 1 ' in the last spin branch (22[. In this 
special case, the spin strings read 

A { q r X=A { q N - 1) +i(N-r + l-2a)c F + 0(i6\c F \), (6) 

with a = 1, . . . , N — r for r = 1, N — 2, respectively. 
In the above equations 5 is a very small number of order 
exp(— L|cf|). 

Substituting the charge bound states k q j with j = 
1,2, ... ,N and the spin strings into Eq. ([3]) results in N 
equations. After multiplying these Af equations together 
and combining with Eq. (U) (see Appendix A) the BAEs 
reduce to 



cx P (MA g L) = TT TT K r f + Z F - (7) 



p=i r=i Aq - A^ - irc F 
The eigenvalues of Hamiltonian ([T]) are then given by 

*2 N N 



E = -N N e b 



2m F 



q=l 



(8) 



with binding energy e b = (h 2 /2m F ) N (N 2 ~ l) c|/12. 

In the strongly attractive limit and in the absence of 
an external field, the A^-fermion clusters are unbreakable 
and we may subtract the binding energy from the energy, 
i.e., 



E F =E + N N e b 



*2 N N 

2m F ^ 9 ' 

q=l 



(9) 



which includes the kinetic energy of the bound clusters 
and the interaction energy produced from cluster-cluster 
scattering. 

In the thermodynamic limit, N F — > oo and L — > oo 
at fixed density n F , the energy of the system can be 
represented in the integral form 



En 



h 2 L 



, A^ / A 2 PF (A) dA 
2m F J_ B 



(10) 



where p F (A) is the density distribution for A determined 
by the integral form of BAEs (|7|) as 



PF (A) 



A^ 



1 N-l rB 



27T TT ' 

r=l 



r \c F \ 



b r 2 c 2 F + (A - A') 



2 p F (A')dA'. 

The integration limit B is determined by the linear den- 
sity n F = N J B B p F (A) dA. 

In terms of the dimensionless energy ejv (l F ), we have 



Ej. 



h 2 L 

2m n 



n%e N {■jp) , 



(12) 
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where — np/N, win = Nmp and 

eN(lF) = ^ / zg N {z)dz. 



(13) 



Here we have defined z = A/B, A = \cf\ /B and gjv (2) = 
p(Bz) /TV. The scaled density distribution g^ (z) is de- 
termined by 



9N (z) 



1 1 

2^ 



1 JV_1 r 1 

r—1 



rA 



5a? (z')dz', 



-i r 2 A 2 + (z — z') 2 
A = N 2 \ 1F \ J g N (z)dz, 

which come from Eq. (|11[) and the linear density rip. 
In the strongly attractive limit with \jp\ ^> 1, we can 
expand the dimensionless energy in terms of 1/| 7 f|. Up 
to 3rd order, this gives 



E F = 



H 2 N N 7T 2 

2m^vi 2 3 
32 fi 



4 12 

iTivl 7/v 



l7Jv| 



2 

15C 3 



(14) 



in which = TV 2 7i?/£, with £ = X^] 1 V r an0 - 



E P =i l n 



A. Equivalence to a super TG gas for TV even 

The unbreakable TV-fermion cluster state is effectively 
described as a composite boson for even TV. Before con- 
structing the mapping relation, we first give a brief re- 
view of the STG state of the attractive Bose gas. The ID 
interacting Bose gas composed of N B bosons with mass 
mg is described by the Hamiltonian 



^ h 2 d 2 



9B 



^T5(xi-Xj), (15) 



2m B dx 2 

2 — 1 1 2<J 



with interaction g B = -2?i /(roBflfn) where is the 
ID s-wave scattering length. The energy eigenvalues are 
given by 



-2 N E 



E F 



2m j 



(16) 



i=i 



where the kj are determined by the BAE [tJ 



exp (ikjL) = — TJ 



kj — ki + ice 
J: kj -fa - ic B ' 

with c B = m B g B /h 2 = -2/ a lD- 



(17) 



For attractive bosons, the ground state solution for 
the BAE (fT7| is a complex st ring solution correspond- 
ing to McGuire's cluster state [ijj- On the other hand, 
the BAE (JTTJ) has real solutions even for c B < 0, which 
correspond to some highly excited states of the attrac- 
tive Bose gas. The super TG state is the lowest gas-like 
state with real solutions for BAE ([17} tULESIHl- In the 
strongly attractive limit, the energy of the STG state of 
the attractive Bose gas can be expressed as 



E 



STG 



h 2 N% 7T ; 

2m s L 2 3 
32 



4 
\ib\ 



12 

\ib\ 



+ - 



\ib\ c 



1 - 
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(18) 



with 7s = c B /n B . 

Comparing equations (TT41 and ([15]) , it is clear that the 
two expressions are identical up to the second order of 
7F if lB = 7JV = A 2 7 F /C, TVg — — Np/N and 
m B = m m = Nwif. Since the TV-bound state formed 
by TV fermions with opposite spin for even TV has a mass 
m B = Nmp, we can conclude that the TV/v TV-bound 
states are equivalently described by the super-TG phase 
of the interacting Bose gas with the effective ID scatter- 
ing length 



1 1D 



C F 

TV 



(19) 



Schematically, we illustrate such a mapping in Fig. 
1(a) by taking the SU(4) Fermi gas as an example. We 
also compare the ground state energy of the SU(TV) 
(TV = 2,4) Fermi gas with the energy of the STG phase 
of the Bose gas composed of composite bosons with mass 
Nitif in Fig. [2J The ground state energy of the SU(2) 
Fermi gas is identical to the energy of the STG phase of 
the Bose gas for all 7^ [Hj], whereas the ground state en- 
ergy of the SU(4) Fermi gas matches very well to that of 
the STG gas for large 7at. As shown in Fig. 3, the rela- 
tive error for the ground state energy of the SU(TV) Fermi 
gas and the energy of the corresponding STG gas is less 
than 0.1 for |7jv| = 10, and less than 10~ 7 for \ \ = 600. 
This indicates that the STG gas provides a good effective 
description for the ground state of the strongly attractive 
SU(TV) Fermi gas, although the mapping is not exact for 
all 7jv like for the SU(2) case. 



B. Equivalence to a super Fermi TG gas for N odd 

The Hamiltonian for the ID p-wave interacting polar- 
ized Fermi gas reads 

Np n 2 d 2 

h p = Y, -^r^2 + 9p E v & ~ ( 2 °) 



2m p dx 2 

where V(x l - x s ) = (gf- - g§-J S( Xi - x s ) (^§- 
is the pseudo-potential for p-wave interaction [25M27 



d 
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FIG. 1: (Color online) The strongly attractive iV-bound state 
Fermi gas can be effectively described by a super Tonks- 
Girardeau gas composed of attractive bosons for (a) even 
N and can be effectively described by a super Fermi Tonks- 
Girardeau gas composed of p-wave repulsive fermions for (b) 
odd N. 



FIG. 3: (Color online) The relative error vs |7jv|, where 
Ae (7jv) = e (7^) — cstg {in)- 

where the quasi-momenta kj are determined by the BAE 
28, 2"' 




FIG. 2: (Color online) Comparison of the ground state ener- 
gies of the attractive SU(iV) Fermi gas and the energy of the 
corresponding STG phase of bosons (for even iV) or p-wave 
fermions (for odd iV) (effective repulsive p-wave fermions) for 
different N with large interaction jn- The inset shows a mag- 
nified view of ground state energies. 



and g p — —2fi a^ D /m p [25l. |26|. The dimensionless in- 
teraction parameter is defined by r ) p = m p g p n p /h 2 = 
~2a^ u n p . For p-wave interacting fermions, the energy 
eigenvalues are given by 



En 



2m t 



N p 



(21) 



exp (ikjL) = — J j 



1=1 



1C K 



(22) 



where the parameter c p = — l/(2a^ D ). 

It is clear that BAE ^2\js identical to BAE dTTJ) 
if c p = cb and N p = Nr |28l. |29| . which means that 
there is a one-to-one correspondence between the p-wave 
Fermi gas and the interacting Bose gas [23| • Correspond- 
ingly, the STG state of the attractive Bose gas has a 
Fermi correspondence which is the lowest gas-like state 
of the weakly interacting p-wave fermions with g p — > + . 
For weakly interacting p-wave fermions in the thermody- 
namic limit, the energy of the lowest gas- like state has 
the form 



2m p L 2 3 



l + 4| 7p | + 12| 7p | 



-:«|l--)| 7p f + 



(23) 



where | 7p | <C 1. 

Comparing equations (fill and (|2"51) , it is clear that the 
two expressions are identical up to the second order in 
j F if 7p = l/ 7Ar = CI (N 2 1F ): N p = N N = N F /N and 
m p = tojv = Nmp. Since the iV-bound state formed by 
N fermions with different hypcrfine states for odd N is 
a composite fermion with a mass m p = NrriF, we can 
conclude that the Nn TV-bound states are equivalently 
described by the gas-like state of the weakly interacting 
p-wave Fermi gas with the effective ID scattering length 



'ID 



AN ' 



(24) 



The mapping is schematically displayed in Fig. 1(b). 
The comparison of the energies of the SU(7V) (N — 3, 5) 
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Fermi gas and the STG phase of the p-wave Fermi gas is 
also given in Fig. [21 which indicates a good matching in 
the limit of large \ jn\- Similarly, as shown in Fig. 3, the 
relative error for the ground state energy of the SU(iV) 
Fermi gas and the energy of the corresponding STG state 
of the p-wave Fermi gas is less than 10 -7 for \jn\ = 600. 
This indicates that, although not exact for all 7jv, the 
mapping provides a very good description for large |7at|. 



III. SUMMARY 

In summary, we have examined the equivalence be- 
tween the ground state of the strongly attractive N- 
component Fermi gas and the super Tonks-Girardeau 
phase of an effective Bose or p-wave Fermi gas. By 
comparing the ground state energy of strongly attractive 
fermions with the energy of the super Tonks-Girardeau 
phase of the Bose gas or p-wave Fermi gas, we find that 
the bound 7V-fermion clusters formed in the strongly at- 
tractive regime should be described by the super Tonks- 
Girardeau phase of attractive composite bosons (for even 
N) or composite fermions with effective p-wave interac- 
tions (for odd N). The super Tonks-Girardeau gas phase 
thus provides an effective description for the ground state 
of strongly attractive multi-component fermions. 
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Appendix A: Derivation of Bethe ansatz equations 
for Fermi bound states 



Here we show how to derive BAEs for Fermi bound 
states by taking as example the system of 3-component 
fermions. For simplicity, we define the function 



e„ (x) 



x + inc f 



x — mcz 



For 3-component fermions the BAEs (J3)) and §4$ are 



Ah 



nei(fc,-AW), (Al) 

Mi 



\8 



p=i 

M 2 



Ile-i AW-Ap) (A2) 



Mi 



1=1 

Mi 



n^AP-AW) = -n*(A, (9, -A£>) (A3) 

a — 1 m=l 

for j — 1,2, ...,Nf, ot,j3= 1,2,..., Mi and l,m = 
1,2,..., M2. Here we confine our attention to the equally 
populated case N F = 3iV 3 , M x = 2N 3 and M 2 = N 3 . 



For strong attraction, i.e., for L\cf\ 3> 1, the charge 
bound states and spin strings are of the form 

k qM = A q +i(4-2hi)c' F + 0(i5\c F \), (A4) 
Ag 2 = Aq+i(3-2h 2 )c' F + 0(iS'\c F \), (A5) 



(21 

for exponentially small S and 5', with A q = Aq , h\ = 
1,2,3, h 2 = 1,2, and q = 1,2,..., A^ 3 . 

For the attractive regime, the common real parts in 
the bound states (|A4[) and spin strings (|A5[) lead to zero 
factors in the BAEs (IAHIA3p . In order to avoid ill-defined 
equations, we eliminate such zero factors in the BAE level 
by level. The first step is to deal with the charge bound 
state kqj^ with hi = 1,2,3 in the BAE ([AT]) , i.e. 



g i3A 9 L _ e i(fc, i i + fc g _ 2 +fcg,3)i 
2 3 



n u^m-^i 



h 2 =l hi=l 
Mi/2 

x Yl e 2 (A g -A Q )e 4 (A (? -A Q ) (A6) 

a=l 



The first terms on the r.h.s. of (|A6[) contain zero factors 
which have to be eliminated. From the second BAE (|A2I) , 
we have 



2 3 



n n«-x(*^-A£ 

h 2 = l h 1 = l 
N F /3 

x Y[ e 2 (A 9 -A Q )e 4 (A 9 -A Q ) 

a=l 
2 



n-x(Aa 2 -A ( 



h 2 = l 
Mi/ 2 



J\ e 2 (A q ~ A Q ) e 4 (A, - A a ) . (A7) 



a=l 



In order to eliminate the first factor on the r.h.s. of 
these equations, we extract this factor from the third 
BAE (fM}, i.e. 

2 Mi/2 

n e -i( A X 2 -A ? ) n e 2 (A g -A a ) 

/i 2 = l a=l 

M 2 

= \{e 2 (At - A m ) . (A8) 



m— 1 



Substituting (KEl into (TAT)) . Eq. becomes 

Mi/2 



A3A„L 



H e 2 (A q - A a ) e 4 (A, - A a ) , (A9) 



a=l 



which is the 3-component fermion form of Eq. ([7]) . 
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